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It has been proved that the charged stringy black holes are stable under the perturbations of 
massive charged scalar fields. However, superradiant instability can be generated by adding the 
mirror-like boundary condition to the composed system of charged stringy black hole and scalar 
field. The unstable boxed quasinormal modes have been calculated by using both analytical and 
numerical method. In this paper, we further provide a time domain analysis by performing a 
long time evolution of charged scalar field configuration in the background of the charged stringy 
black hole with the mirror-like boundary condition imposed. We have used the ingoing Eddington- 
Finkelstein coordinates to derive the evolution equation, and adopted Pseudo-spectral method and 
the forth-order Runge-Kutta method to evolve the scalar field with the initial Gaussian wave packet. 

It is shown by our numerical scheme that Fourier transforming the evolution data coincides well with 
the unstable modes computed from frequency domain analysis. The existence of the rapid growth 
mode makes the charged stringy black hole a good test ground to study the nonlinear development 
of superradiant instability. 
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I. INTRODUCTION 

Superradiance is an interesting classical effect of the 
rotating or charged black holes [lj- 4]. Although in any 
classical process in black hole spacetime no particle (even 
for photon) can escape from the event horizon of black 
hole, the energy of rotating or charged black hole can be 
extracted out via the superradiance process. It should 
also be noted that this classical superradiance process 
does not violate the area theorem. Consider, for example, 
the superradiance effect of rotating black hole. Suppose 
a bosonic wave of the form with u> and m 

being the frequency and azimuthal quantum number, is 
scattered by the event horizon. If Re[u>\ < mUln with fin 
being the angular velocity of event horizon, the incident 
wave gets amplified by the black hole rotating. For the 
charged black holes, the charged bosonic fields should be 
considered in order to create this effect, where the charge 
of perturbation field play the role of rotating [g] . 

This effect may trigger an instability if the wave re- 
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peatedly scatters off the black hole, which is named as 
superradiant instability or black hole bomb mechanism 
[6, 7|. This can be achieved by either introducing a mass 
term of the bosonic field or adding a mirror-like bound¬ 
ary condition outside the black hole, which will reflect 
the bosonic wave repeatedly. Then the energy extracted 
from the black hole can grow exponentially with the time. 
For many cases, this instability triggered by superradi¬ 
ance has been confirmed by both analytical and numeri¬ 
cal methods [§-|33|]. One can refer to Ref. 0] for a recent 
review on this topic. However, this is not the full story. 
It is shown that, in some cases, the scalar fields do be¬ 
have the superradiance effect, while the instability does 
not occur. For example, the Reissner-Nordstrom black 
hole is stable against the sunerradiance [351438! . 

In [ 39 ), we have shown that the mass term of charged 
scalar field can not support quasi-bound state outside 
the charged stringy black hole. This is to say that the 
charged stringy black hole is stable against the massive 
charged scalar perturbation. This result was further con¬ 
firmed in pQt by computing the quasinormal modes of 
charged scalar perturbation. In [4l[ , we have studied the 
superradiant instability of the scalar field in the back¬ 
ground of charged stringy black hole due to a mirror¬ 
like boundary condition. The analytical expression of 
the unstable superradiant modes is derived by using the 
asymptotic matching method Soon later, this analyt- 
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ical result is further confirmed by the recent numerical 
studies 0 . Motivated by the work 0 , the station¬ 
ary cloud-like configurations of scalar field satisfying the 
superradiant critical frequencies are also studied in 
However, the instability is only revealed by computing 
the unstable boxed quasinormal modes in the frequency 
domain analysis (55j . 

In this paper, we shall perform the time domain anal¬ 
ysis of the charged scalar field configuration in the above 
charged stringy black hole background. Although the 
time domain analysis of superradiant instability is a chal¬ 
lenging work due to the small growth rate, a time do¬ 
main study of superradiant instability for the scalar field 
on Kerr spacetime has been initiated in |45j]. Instead of 
the usually used Schwarzschild like coordinates, we will 
exploit the ingoing Eddington-Finkelstein coordinates, 
where the partial differential evolution equation is gener- 
ically first order in time such that there is no need for one 
to introduce an auxiliary field for numerical computation. 
In addition, compared to the finite difference method, the 
pseudo-spectral method has become an increasingly pop¬ 
ular numerical technique to study the black hole physics 
in particular due to its higher computation efficiency. For 
example, the pseudo-spectral method has recently been 
used to calculate the quasinornral modes of asymptoti¬ 
cally flat black holes in [46] and to study the gravita¬ 
tional wave collapse in 0 |. One can refer to Ref.jHj 
for a review of this numerical method. With this in 
mind, we shall resort to the pseudo-spectral method in 
the spacial direction supplemented with the forth order 
Runge-Kutta method in time to numerically evolve an 
initial Gaussian wave packet. As a result, the unstable 
modes, obtained from the Fourier transforming the evo¬ 
lution data, prove to be in good agreement with those 
computed from the frequency domain analysis. We also 
confirm the existence of the rapid growth modes, which 
makes the charged stringy black hole a good test ground 
to study the nonlinear development of superradiant in¬ 
stability. 

This paper is organized as follows. In Sec. II, we 
describe the basic setup of a charged scalar field with 
mirror-like boundary condition in charged stringy black 
hole background. The involved partial differential equa¬ 
tion that governs the evolution is also derived. In Sec. 
Ill, our numerical scheme is detailed. We present our nu¬ 
merical results and relevant discussions in Sec. IV. Final 
remarks are made in the last Section. 


II. FRAMEWORK 


The black hole background considered in the present 
paper is a static spherically symmetric charged black hole 
in the low energy effective theory of heterotic string the¬ 
ory in four dimensions, which was firstly found by Gib¬ 
bons and Maeda in [0 as well as independently found by 
Garfinkle, Horowitz, and Stronringer in p50j a few years 


later. The metric is explicitly given by 


ds 2 



dt 2 + 



-l 


dr 2 


( dd 2 + sin 2 9d(j) 2 ) , 


(1) 


and the electric potential and dilaton field are given by 



2<s> = l _ Q_ 
Mr 


( 2 ) 


The parameters M and Q are the mass and electric 
charge of the black hole, respectively. The event hori¬ 
zon of black hole is located at r = 2 M and the singu¬ 
larity is located at r = Q 2 /M. We only consider the 
non-extremal case in which the parameters satisfy the 
condition Q < \f2M. 

For the latter convenience of numerical calculation, we 
would like to adopt the ingoing null coordinate defined 

by 


v = t + r* , 


( 3 ) 


with r* being the tortoise coordinate 


* 


r 


dr 

2 M 


= r + 2M In 


2 M 1 


( 4 ) 


Then the metric in the ingoing Eddington-Finkelstein co¬ 
ordinates reads 


ds 2 = — ^1 — dv 2 + 2 dvdr 

+r ^ 2 S ™ 2 ’ (5) 

and the electric potential can be written as 

A v = -Q, ( 6 ) 

r 

where we have also performed a gauge transformation 
A-Vij) with ip = Q f 

For simiplicity, it is sufficient to consider the test 
charged massless scalar field, the dynamics of which is 
governed by the Klein-Gordon equation 


{V u -iqA v ){V v -iqA v )^ = 0 , (7) 


where q denotes the charge of the scalar field. 
By taking the ansatz of the scalar field 


l,m 


l 

s/r(r — 2 M) 


<tnm{u,r)Y lrn (6,(j)) , 


( 8 ) 


the Klein-Gordon equation is shown to be separable in 
this background. After some algebra, we can finally get 
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the following partial differential equation 

Q 

Q 4 r — 2 M 


, r ' 2M a2i *«Q a i if lQ 

d v dr<p H-——+ ^-<+-0 + —^dr<P - -^J<; 


1 

'2 r 3 (r~Q 2 /M) 

1(1 + 1 ) , 

2 r(r-Q 2 /My 


Q 2 - 2 Mr + 


= 0 , 


4 M 2 r-Q 2 /M 


(9) 


where we have dropped the mode number subscripts to 
simplify the notation. 

We can now evolve the scalar field according to this 
partial differential equation numerically with the initial 
condition. From the frequency domain analysis, we know 
that the mirror-like boundary condition should be im¬ 
posed to generate the superradiant instability of charged 
scalar filed in the stringy charged black hole spacetime. 
This condition now can be explicitly expressed as 


4>(v,r m ) = 0 , 


( 10 ) 


where r m is the location of mirror. This condition should 
be satisfied during the whole time evolution process. 


III. NUMERICS 


We use the method of lines in the forth order Runge- 
Kutta scheme to integrate in time. The space direction is 
then approximated by using the pseudo-spectral method 
instead of the finite difference method. To be more spe¬ 
cific, we first use a linear transformation 


2 (r - 2 M) 
r m - 2 M 


( 11 ) 


to map the r coordinate interval [2 M, r m \ into an inter¬ 
val [—1,1]. We then expand the scalar field in a basis 
of Chebyshev polynomials in the space direction. The 
scalar field function is discretized by the values of the 
expansion on the Gauss-Lobatto collocation points. The 
spatial derivative of the scalar field is computed by multi¬ 
plying the Gauss-Lobatto derivative matrix. In the time 
direction, we use the standard fourth order Runge-Kutta 
scheme to evolve the scalar field. To be more precise, 
we first evolve the space derivative of the scalar field ac¬ 
cording to HI and then integrate out the scalar field by 
multiplying the integral matrix with the built in mirror¬ 
like boundary condition (HI. Therefore the mirror-like 
boundary condition is automatically maintained in the 
whole evolution process. The above procedure can be 
realized in the software of MATHEMATICA. It is note¬ 
worthy that the similar numerical trick has also been used 
in [5l| to study the holographic vortex dynamics in the 
context of AdS/CFT correspondence. 


IV. RESULT 

Figure 1 shows a generic evolution of scalar field. The 
initial configuration at v = 0 is a Gaussian wave packet 


TABLE I: Boxed quasinormal modes for the first 2 mirrored 
states with M = l,r m = 30,1 = 1,Q = 0.6, q = 0.6. The 
second column is taken from the frequency domain analysis, 
while the third column is read off from the Fourier transfor¬ 
mation of evolution data, i.e. the third panel in Figure 1. 


n 

UJ 

Re[uj\ 

0 

0.158905 + 5.69687 x 10 ~'i 

0.159 

1 

0.254763 - 1.79498 x 10 -5 i 

0.255 


TABLE II: Boxed quasinormal modes for the first 2 mirrored 
states taking M = l,r m = 30, l = 0 ,Q = 0.6, q = 0.6. The 
second column is taken from the frequency domain analysis, 
while the third column is read off from the Fourier transfor¬ 
mation of evolution data, , i.e. the third panel in Figure 2. 


n 

UJ 

Re[cj\ 

0 

0.119962 + 4.85725 x 10 

0.120 

1 

0.211098 - 7.12823 x 10 4 i 

0.212 


centered at r cg , i.e., 

(r-rcg) 2 

0(0 ,r)=(j) 0 e ^ , (12) 

with <j>o = 3 x 10 -5 , r cg = 15 M, a 2 = 30 M 2 . From the 
first panel in Figure 1, we find that at early stage of evo¬ 
lution the amplitude of scalar field oscillatorily decays 
with the time. This is mainly caused by the absorp¬ 
tion of the black hole horizon. While, in the following 
stage, it oscillatorily grows very slowly. The reason is 
that the imaginary part of unstable mode in this case 
is very small, which is shown in Table I. The second 
panel shows the logarithm of scalar field amplitude evolv¬ 
ing with the time. The third panel is the corresponding 
Fourier transformation of the evolution data, where some 
peaks show up at a set of distinct frequencies with the 
unstable mode dominant. As further demonstrated in 
Table I, the magnitude of these distinct frequencies is in 
good agreement with the boxed quasinormal modes in 
the frequency domain analysis. We find the consistency 
between the modes from two sides. The growth time 
scale is set by the inverse of the imaginary part. So only 
after a long time evolution of order of v/M ~ 10 7 can a 
significant growth of scalar field amplitude be observed 
in this case. In what follows we will focus on the case in 
which the unstable mode has a much bigger imaginary 
part. 

By calculating the boxed quasinormal modes, we find 
that the imaginary parts of l = 0 unstable modes are 
generically three orders of magnitude greater than that 
of l = 1 unstable modes. This indicates that the l = 0 
unstable modes have a much faster growth rate than the 
1 = 1 unstable modes, which is confirmed by Figure 2 for 
the real time evolution of the l = 0 case. In particular, 
the exponential growth rate can be estimated as 5 x 1CP 4 
or so by the time evolution of logarithm of scalar field 
amplitude, which coincides with the imaginary part of 
the unstable mode in Table II. Furthermore, as shown in 
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FIG. 1: Time evolution of scalar field with the initial Gaussian wave packet. The parameters are taken as r m = 30, l = 1, Q = 
0.6,and q = 0.6, with M = 1 as our unit. The first panel shows the real part of the scalar field evolving with the time, while 
the second panel shows the logarithm of scalar field evolving with the time. Fourier transformation of the time evolution data 
is demonstrated in the third panel. 



FIG. 2: Time evolution of scalar field with the initial Gaussian wave packet. The parameters are taken as M = 1, r m = 30, l = 
0, Q = 0.6, and q = 0.6. The first panel shows the real part of the scalar field evolving with the time, while the second panel 
shows the logarithm of scalar field evolving with the time. Fourier transformation of the time evolution data is demonstrated 
in the third panel. 


Table II, the real part of unstable mode read off from the 
Fourier transformation of evolution data is also consistent 
with the frequency domain analysis. 

In the above two cases, we have chosen the parame¬ 
ters to make the black hole-nrirror system possess only 
one unstable mode. Note that the superradiant critical 
frequency is given by 152 


The critical frequency for the above two cases, sepa¬ 
rating growing modes (w < cu c ) from decreasing modes 
(w > w c ), is given by 0.18, consistent with the above nu¬ 
merical results. In order to see more unstable modes show 
up, it is reasonable to increase the critical frequency by 
cranking up either the charge of the black hole or scalar 
field. 

We first compare the time evolutions of scalar field 
for different values of black hole charge Q. As discussed 
above, in order to observe the distinct growth of scalar 
field amplitude, we have also chosen the parameter 1 = 0. 
Keeping M = 1 ,r m = 30, and q = 0.6, in the following 
we take three cases with Q = 0.8,1.1, and 1.4 on top of 
the previous case of Q = 0.6. The superradiant critical 
frequencies in these three cases are 0.24, 0.33, and 0.42, 
individually. In Figure 3, we display the time evolutions 
of the scalar field amplitude (left three panels). The first 
3 boxed quasinormal modes are displayed in Table III. As 
expected before, the number of unstable modes increases 
with the black hole charge. When Q = 0.8 and Q = 1.1, 
there are two unstable modes. When Q = 1.4, there 
are three unstable modes. These unstable modes are all 
displayed in the Fourier transformation of the evolving 


data (right three panels). It is also shown that the mode 
dominates the evolving is always the one with the largest 
positive imaginary part in Table III. However, there is no 
general rule of the instability of black hole-mirror system 
depending on the black hole charge Q. In other words, 
the fast or the slow growth of scalar field is not deter¬ 
mined by the large or the small black hole charge. 

Next, we compare the time evolutions of scalar field for 
different values of the field charge q , keeping the same ini¬ 
tial Gaussian wave packet as before with M = 1 ,r m = 
30,1 = 0, and Q = 0.8. In Figure 4, we display the 
time evolutions of the field amplitude (left three pan¬ 
els) and the Fourier transformation of the evolution data 
(right three panels). The boxed quasinormal modes for 
the first 3 mirrored states are displayed in Table IV. The 
superradiant critical frequencies in these three cases are 
0.16, 0.26, and 0.40, individually. As expected, the num¬ 
ber of unstable modes also increases with the charge of 
scalar field. When q = 0.4, there is only one unstable 
mode. When q = 0.65, there are two unstable modes. 
When q = 1.0, there are three unstable modes. Sim¬ 
ilarly, the Fourier transformation of the evolution data 
shows that the dominant mode is always the one with 
the largest positive imaginary part in Table IV. In ad¬ 
dition, consistent with the previous analytic result pil] ]. 
we also observe that when the field charge q increases, 
the black hole-mirror system becomes more unstable and 
exhibits a more rapidly growing behavior. 

We also show the time evolutions of scalar field for 
different values of mirror radius r m in Figure 5. The cor¬ 
responding boxed quasinormal modes for the first 3 mir¬ 
rored states are also displayed in Table V. As expected, 
Fourier transformation of time evolution data is in good 
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TABLE III: Boxed quasinormal modes for the first 3 mirrored states taking M = 1, l = 0, r m = 30, q = 0.6. 

n Q = 0.8 Q = 1.1 Q = 1.4 

(! 0T28718 — 7". (h 59 ■ in H 04 11096 + 844972 .Jo H 045229 - 6)23948 - 10 W 

1 0.223744 + 2.33795 x 10 _4 i 0.24119 + 5.9991 x 10" 4 * 0.257199 + 1.24213 x 10" 4 i 

2 0.309293 - 2.52535 x 10 _3 i 0.331962 - 2.75849 x 10 _5 i 0.352634 + 4.42053 x 10" 5 i 


TABLE IV: Boxed quasinormal modes for the first 3 mirrored states taking M = 1, r m = 30, l = 0, Q = 0.8. 

n q = 0.4 q = 0.65 q = 1.0 

ii 0ll7342TT09664 > !() H (i. 131152 + 843298 in “ Ol 19771 - L62917 It) *T 

1 0.207633 - 1.04707 x 10 _3 i 0.227552 + 4.26772 x 10 -4 * 0.252281 + 1.29866 x 10 -3 * 

2 0.289823 - 7.09925 x 10 _3 i 0.314069 - 1.76193 x 10 _3 i 0.345103 + 8.61385 x 10 _4 i 



FIG. 3: Time evolutions of scalar field for different black hole 
charges. From top to bottom, Q = 0.8,1.1, and 1.4, individu¬ 
ally. The parameters are taken as M = 1, r m = 30, l = 0, and 
q = 0.6. The left three panels show the real part of the scalar 
field evolving with the time, while the right three panels show 
the Fourier transformation of the time evolution data. In fact, 
the last panel shows the logarithm of Fourier transformation, 
otherwise the other modes would be much overshadowed by 
the most unstable mode. 

agreement with the frequency domain analysis. 

It is shown in the work [53j, (54] that there exists 
rapid growth unstable modes in the charged Reissner- 
Nordstrom black hole. Motivated by this, we would 
like to end this section by demonstrating a very rapidly 
growing time evolution in our model. As such, we take 
the parameters as M = l,r m = 5,/ = 1,Q = 1, and 
q = 30, for which the most unstable mode is given by 
to = 8.01818 + 0.05738571. Figure 6 displays such a time 
evolution of scalar field the Fourier transformation of the 
evolution data for the initial Gaussian wave packet 

. (r —2.5 M) * 1 2 

(f>(0 ,r) = 3 x 10 5 e o.o3M^ . (14) 



FIG. 4: Time evolutions of scalar field for different values of 
the field charges. From top to bottom, q = 0.4, 0.65, 1.0, indi¬ 
vidually. The parameters are taken as M = 1, r m = 30 ,1 = 0, 
and Q = 0.8. The left three panels show the real part of the 
scalar field evolving with the time, while the right three panels 
show the Fourier transformation of the time evolution data. 
The last panel also shows the logarithm of Fourier transfor¬ 
mation. 


Obviously the real part of scalar field evolves with the 
time very rapidly. Figure 6(a) and (c) represent the 
growth of the real part of scalar field and its logarithm. It 
can be easily estimated from Figure 6(c) that the growth 
rate is about 0.05, which is in good agreement with the 
imaginary part of the above unstable mode. In addition, 
we also zoom in on the time evolution to illustrate the os¬ 
cillating behavior of the scalar field in Fig. 6(b). To read 
off the precise value of this oscillating frequency, we fur¬ 
ther take the Fourier transformation of the evolution data 
and the resultant real part of dominant unstable modes 
is 8.015, which is also consistent with the frequency do¬ 
main analysis as well. Our numerical results show that 
the rapidly growing evolution process is dominated by 
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TABLE V: Boxed quasinormal modes for the first 3 mirrored states taking M = 1, l = 1, Q = 1, q = 2. 

Tl I'm — 20 V m — 15 V m —— 10 

() 0365882TT48242 - in H CU7635TX752T6 ■ in 1 i 067915+^08701 - 10 *T 

1 0.53467 + 1.13535 x 10” 4 * 0.686655 + 1.29391 x 10 _4 i 0.947995 + 1.0283 x 10~ 4 * 

2 0.684008 + 9.61229 x 10" 5 i 0.868147 + 1.01999 x 10“ 4 i 1.15524 - 4.12335 x 10“ 3 i 







j 


- 

L _ 

l • 



FIG. 5: Time evolutions for different values of mirror radius. 
From top to bottom, r m = 20,15, and 10, individually. The 
parameters are taken as M = 1,Z = 1,Q = 1, and q = 2. The 
left three panels show the real part of the evolving scalar field, 
while the right three panels show the Fourier transformation 
of the evolution data. 



(a) (b) 



(c) (d) 

FIG. 6: Rapidly growing time evolution. The parameters are 
taken as M = 1, r m = 5, l = 1, Q = 1, and q = 30. (a) Real 
part of the scalar field evolving with the time; (b) Oscillating 
behavior of real part; (c) Logarithm of scalar field evolving 
with the time; (d) Logarithm of Fourier transformation of the 
evolution data. 


the most unstable mode and other unstable modes are 
much suppressed. 


V. CONCLUSION 

In this paper, we have presented a time domain anal¬ 
ysis of long time evolution of charged scalar field con¬ 
figuration in charged stringy black hole background by 
imposing the mirror-like boundary condition. The nu¬ 
merical scheme we have used is pseudo-spectral method 
in the space direction supplemented with Runge-Kutta 
method in the ingoing Eddington-Finkelstein time direc¬ 
tion. As a result, this numerical scheme turns out to be 
highly efficient and amenable to a long time evolution. 
In particular, our numerical results confirm the superra¬ 
diant instability of scalar field in this background when 
the black hole is enclosed in a cavity. Furthermore, it 
is shown that Fourier transforming the evolution data 
coincides with the unstable modes calculated from the 
frequency domain analysis. This paper may provide an 
alternative method to study the time evolution question 
of field perturbations in black hole background. 

In the present work, the backreaction of scalar field on 
the background geometry is not considered. During the 
superradiant process, the field extracts energy from the 
black hole and its energy grows exponentially. Then the 
nonlinear interaction between the field and black hole is 
required to be taken into account in particular for the 
case in which there is a very rapidly growing mode. We 
demonstrate there exists the rapid growth mode exist in 
charged stringy black hole-mirror system. To see what 
final state this superradiant instability leads to, the fully 
nonlinear evolution of the coupled system of scalar field 
and black hole is required, although it is supposed to 
much involved. We hope to address this challenging issue 
in the near future. 
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